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In the recent experiments [Chmiola et al, Science 313, 1760 (2006); Largeot et 
al, J. Am. Chem. Soc. 130, 2730 (2008)] an anomalous increase of the capaci- 
tance with a decrease of the pore size of a carbon-based porous electric double-layer 
capacitor has been observed. We explain this effect by the image forces which ex- 
ponentially screen out the electrostatic interactions of ions in the interior of a pore. 
Packing of ions of the same sign becomes easier and is mainly limited by steric in- 
teractions. We call this state 'superionic' and suggest a simple model to describe it. 
The model reveals a possibility of a voltage-induced first-order transition between a 
cation(anion)-deficient phase and a cation(anion)-rich phase which manifests itself 
in a jump of capacitance as a function of voltage. 

PACS numbers: 82.47.Uv 



* Electronic address: skondrat@ic.ac.uk 

t Electronic address: a.kornyshev@imperial.ac.uk 



2 



Electric double layer supercapacitors (EDLCs) have recently attracted a considerable 
attention due to their high energy density and fast power delivery jl, 2|. The carbon-based 
EDLCs store energy at the electrolyte/carbon interface, and in order to increase the amount 
of the energy stored it is necessary to increase the carbon surface area. This is usually 
achieved by using highly porosous carbon materials, with the pore size being twice the size 



of the ions to allow ions adsorp 
In the recent experiments 



;ion on both pore walls 



aa. 



y, |5| Gogotsi and co-workers have measured the capacitance 
per surface area of EDLCs with nanoporous carbon-based electrodes vs the pore width. They 
obtained a surprising result: an anomalous increase of the capacitance with decreasing the 
pore width for nano and subnano pores comparable in size with the size of bare ions. Both 
an organic electrolyte J4] and a (solvent-free) ionic liquid [5] have been tested as electrolyte 
medium leading to the same effect. 

Molecular Dynamics (MD) simulations have confirmed the experimental results 
For instance, in Ref. [3] the simulations have been performed for EMI cations and TSFI 
anions as ionic liquid and the interior of a single-wall carbon nanotube (CNT) as pore. The 
'anomalous' increase of the capacitance has been observed for the CNT radii down to the 
ions size. 

It has been proposed that partial desolvation of the ions entering small pores is responsible 
for the anomalous behaviour of the capacitance in the case of organic electrolytes (4]. A 
recently considered electric wire-in-cylinder model 8J offers a good fit to the experimental 
data but for further development of the porous EDLCs it is important to understand the 
underlying physics of the capacitance in nanopores. 

In this Letter we propose a simple explanation of the anomalous capacitance behaviour by 
introducing the concept of a superionic state in a metallic nanopore caused by (i) exponential 
screening of the ion-ion pair interactions in a nanogap, which arises due to the image forces 
and enables more ions of the same sort to occupy the pore. By the same token, there is 
also (ii) an electrostatic contribution to the free energy of transfer of ions from the bulk to 
the pore interior. The latter tends to increase the ion density in a pore, compensating the 
opposing trend due to the loss of solvation shell by the ions moving into the pore. 

Let us consider a point charge at z — z\ confined between two metal plates (c.f. Fig. [1] 
with V — 0). The electrostatic potential satisfying the Laplace equation and Dirichlet 
boundary conditions at the plates can be easily found using the Fourier-Bessel transform to 
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FIG. 1. Schematic view of the cross section of a single, laterally infinite, slit-like narrow pore as 
a part of a porous electrode. The pore width is L, V is the voltage with respect to the reference 
electrode (in the "bulk" outside the pore), and d is the diameter of ions. 
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where e is the dielectric constant of a medium between the plates, R is the lateral distance 
from the charge, L is the plates separation (i.e., the pore size), and K n (x) is the modified 



Bessel 
x > 1 



'unction of the second kind of order n. Using the asymptotic behaviour of K Q (x) for 
9j one finds for R 3> L/tt 
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sm(irzi/L) sm(7rz/L), 



(2) 



eV2LR 

which means that the electric potential exponentially decays with the distance in the plane 
parallel to the plates, and the Coulomb interaction between two charges is effectively screened 
out. In the opposite limit of large L ^> irR, Eq. ([T]) recovers the Coulomb law. 

We consider for simplicity a single, metallic [lOj, slit-like 11] pore, infinitely extended in 
the lateral directions (see Fig. [1]), and neglect the boundary effects at the "entrance" and 
"exit" or closing of the pore (for supercapacitor electrodes the typical pore length is usually 
much longer then its width). The electrostatic potential at the pore walls is measured 
with respect to the reference electrode in the 'bulk' of the capacitor. We consider a small 
nanoscale pore, in which there is room only for one ionic layer (as in the experiments and 
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simulations). Furthermore, we assume for simplicity that the ions are situated on average 
in the middle of the pore (i.e., z — L/2, c.f. Eqs. (J3J) and (TjJ). Although not generally true, 
this assumption should not qualitatively change our conclusions. 

Using now Eq. (pQ), and taking the voltage into consideration, one finds for the internal 
energy (per surface area) due to electrostatic interactions of an ionic liquid in a pore 

(3U(p,c) = f3eVc + 4c 2 R c (p)L B J2 Sm ™ ] - K x {vmR c {p) / 'L) , (3) 

m=l 

where /3 = l/k B T, with k B being the Boltzmann constant and T the temperature, and 



Lb = (3e 2 /e is the Bjerrum length [12J; c = Z + p + — Z_p_ and p = p + + p_ are the two- 
dimensional charge density and total density of the ions in the pore, respectively, where p± 
are the (two-dimensional) densities of the ± ions and Z± the corresponding valencies. In 



the derivation of Eq. we have used the cut-out disk approximation [13j with the cut-out 
radius R c (p) = (vrp)^ 1 / 2 . The voltage, which appears linearly in Eq. ([3]), drives the anions 
to and cations from the pore for V > 0, and vice versa for V < 0. 

We present the free energy of transfer of an ion from the bulk into a pore as a sum of 
two contributions: (i) a "re- solvation" free energy 5E a (a = {+,—}) due to the (partial) 
loss of the solvation shell when an ion moves from the bulk to the pore, which we assume to 
be L-independent, and (ii) a change in the self-energy of an ion due to confinement, which 
characterizes its interaction with the metal plates. The latter can be approximated as the 
work of "charging" of a point charge inside the pore minus the work needed to charge it in 
the bulk with the same dielectric constant: (q 2 /2) lim r _ 5 . ri (0 i (r, r x ) — 4>l=oo{ v , r i))- Using 
Eq. (pQ) one obtains for this contribution to the internal energy (per surface area) 

0E S ( P± ) = (PS E * - l b/L /(1/2) zi) Pai (4) 

a=± 

where 

/w-fG- ^y )"' 

The function f(x) is positive-definite around f(x = 1/2) = min xg [ 0;1 ] f{x) = ln(2), and 
hence the second term in Eq. (j2J) is negative. The resolvation energy 5E± is positive and 
depends on the type of electrolyte medium; its value is not known exactly but we estimate 
it to be not larger than 20 k B T jjjjj.This implies that for sufficiently large L B /L the total 
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free energy of transfer is negative, and will therefore favour an increase of the ion densities 
inside the pore relative to the bulk (irrespective of the voltage and the sign of the charges). 
The total free energy of an ionic liquid in the pore (per surface area) is 

F=U + E a + ^2fji a p a -TS, (6) 

a=± 

where S is the entropy and p± chemical potentials. For the ions of the same size d the 
entropy is given by 15] 



where Vq = 7id 3 /Qr] m is the minimal volume per ion and rj m the maximum packing frac- 
tion [16j . The first term in Eq. (j7J) corresponds to the entropy of ions, and the last term is 
the entropy of a solvent (if there is one). Note that Eq. (JZJ) reflects the three-dimensionality 
of the problem [l7J. Finally, the chemical potentials /i± are set to the chemical potentials 
of the ionic liquid in the "bulk", i.e., where the electric potential vanishes (we recall that, 
besides here, the bulk part does not appear in the model); one gets 

-(oo) \ / rv 

v oP± \ , 1 Z± 



where Z — Z + + Poo = p^ + p^ is the total (three-dimensional) density of the ions in 



the bulk, 7 = foPoo (see Ref. Il5|), and we have used the electro- neutrality condition obeyed 
in the bulk. 

The free energy fl6]) can now be minimized numerically with respect to the ion densities 
p±, and the differential capacitance can be found by numerically differentiating the negative 
of the total charge accumulated in the pore with respect to the voltage. 

The capacitance (per surface area) vs the pore width is shown in Fig. [2]^a) for a few 
values of the voltage V. At zero voltage the capacitance increases with decreasing the pore 
width (solid line in Fig. I2fa)), in agreement with the experimental observations [4, |5|. Note 
that C > Ch for small pore widths, where Cu = e/2ird is the Helmholtz capacitance 18 |. 
Obviously, C = for pores smaller than the ion size because the ions cannot enter the pore 
(not shown in Fig. [2]). The differential capacitance also vanishes for voltages higher than a 
certain threshold voltage V (L) (see Fig. [2(b) and dot-dash line in Fig. |2(a)); this is because 
the pore is fully occupied by anions for V > V . 
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FIG. 2. (color online) Differential capacitance per surface area, expressed in terms of the Helmholtz 
capacitance Ch = s/2nd, as a function of (a) pore width L/d and (b) voltage V/Vr, where 
Vt = l//3e is the thermal voltage (~ 26 mV at room temperature). The ion diameter is d = 0.7 nm, 
Bjerrum length Lb = 20 nm, r] m = tt/6, 7 = 0.5, re-solvation energies 5E± = WksT, and valencies 
Z± = 1. The capacitance increases with decreasing the pore width at zero voltage (solid line in (a)), 
in accordance with the experimental observations (see Refs. P, 5]). There is a jump in capacitance 
at L/d m 1.63 for V/V T = 20 in (a) and at V/V T ~ 25.2 for L/d = 1.2 in (b) manifesting a first 
order transition. A maximum for L/d = 1 in (b) reflects a close vicinity of a critical end-point 
(c.f. Fig. El). 



For larger pores and intermediate voltages we observe a jump in both capacitance vs pore 
size and capacitance vs voltage curves (dash lines in Fig. [2]). It marks the onset of a first 
order phase transition. This transition occurs when the applied voltage (positive in this 
case) is sufficiently high to drive the cations out of the pore but is not yet high enough to fill 
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FIG. 3. (color online) (a) The total density p and the cation density p+, expressed in terms of 
Pmax = L/vq, as a function of voltage for the pore width corresponding to a first order transition 
(L/d = 1.2) and a critical end-point (L/d ~ 1.08). (b) The phase diagram spanned in the (voltage, 
pore size) plane. The solid line corresponds to a first order transition between a cation-deficient 
phase and a denser, cation-rich phase, and the full circle denotes a critical end-point. Thin vertical 
lines correspond to the profiles in (a). The model parameters are the same as in Fig. [2j 



the pore completely with anions. In other words, the transition separates a cation-deficient 
(CD) phase and a denser cation-rich (CR) phase. (In the case V < the transition separates 
the anion-rich (AR) and anion-deficient (AD) phases.) This is demonstrated in Fig. |3ja) 
where we plot the cation density p + vs voltage. At L/d = 1.2 the density profile p+{V) 
is discontinuous, manifesting a first order transition. At smaller L/d ~ 1.08 the profile is 
continuous but dp + /dV diverges at V/Vr ~ 24.8 which indicates a critical end-point. The 
corresponding phase diagram in the (voltage, pore-size) plane is shown in Fig. [31(b) . 
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The CR/AR phases are characterized by high total ion density close to p max = L/v (see 
Fig. 13(a)). This is due to the superionic state created in the pore interior, i.e., the favourable 
free energy of transfer [see Eq. (j3J)] and the exponential screening of the Coulomb interactions 
[see Eqs. fl2]) and (J3J)]. The pore charging in this case is mainly due to the exchange of the 
cations in the pore with the anions from the outside of the pore (for V > 0) such that the 
total density remains practically constant (see Fig. 13(a))- This is to be contrasted with the 
charging in the CD and AD phases, which is essentially characterized by an increase of the 
total ion density. Therefore, the phase transition lines also separate two different 'charging 
regimes.' It is reasonable to expect the charging/ discharging dynamics to be also different in 
these two regimes, which may be of a considerable importance for power delivery of porous 
EDLCs. 

In summary, we have considered a simple phenomenological model in order to under- 
stand the properties of EDLCs with nanoporous electrodes. Our model takes into account 
the exponential screening of the ion-ion interactions in a metallic nanopore and the inter- 
action of ions with pore walls - both determined by image forces, which underpin what we 
call a 'superionic state' in metallic nanopores. It explains an 'anomalous' increase of the 



capacitance with a decrease of the pore size {4, 5|. We have also predicted an interesting 
voltage-induced first-order transition between a cation(anion)-deficient phase and a denser 
cation(anion)-rich phase, manifested as a jump in the capacitance. In spite of the coarse 
nature of the model and a number of simplifying assumptions, we believe it leads to qual- 
itatively correct results. It would be beneficial to check its predictions using, for instance, 
a more robust density-functional theory or Monte Carlo/Molecular Dynamics simulations, 
which can also probe the discussed effects in a wider range of pore sizes, and of course to 
verify them by experiments. 
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